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Abstract. The multiplicity function (MF) of groups and clusters of galaxies is determined using galaxy catalogues 
extracted from a set of Digitized Palomar Sky Survey (DPOSS) plates. The two different types of structures (of low 
and high richness) were identified using two different algorithms: a modified version of the van Albada method 
for groups, and a peak finding algorithm for larger structures. In a 300deg 2 area up to z < 0.2, we find 2944 
groups and 179 clusters. Our MF covers a wide range of richnesses, from 2 to 200, and the two MF's derived by 
the two algorithms match smoothly without the need for additional conditions or normalisations. The resulting 
multiplicity function, of slope a — —2.08 ± 0.07, strongly resembles a power law. 
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1. Introduction 

The multiplicity function (hereafter MF), in its differen- 
tial form, is defined as the number of groups or clusters 
per area or volume unit and per richness unit. 
The MF, which is the richness spectrum of galaxy aggre- 
gates, parametrises the observed clustering of galaxies and 
hence, together with the correlation and luminosity func- 
tions, is one of the fundamental cosmological observables. 
With respect to the complete description of clustering, the 
MF is complementary to the covariance function (which 
is related to the two-point correlation function) , being re- 
lated to the ratio of the amplitude of the higher-order to 
the two-point correlation functions l|Gott Ez Turneilll977l 
hereafter GT). Due to computational costs and errors, the 
measurement of correlation functions of order N becomes 
unreliable for N > 3, and the MF is therefore a crucial 
means of obtaining information on higher order cluster- 
ing 



The Press-Schechter theory l|Press fc Schecterl Il974|) 
states that the shape of the mass function (a power-law 
mass distribution with an exponential cutoff at the bright 
end) should provide important clues concerning the condi- 
tions at the epoch of recombination and does not depend 
on the cosmic density parameter f2. The steepness of the 
initial density fluctuation spectrum constrains the broad- 
ness of the mass function. 
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The MF, the mass function or the luminosity function all 
describe in a similar way the cosmic abundance of objects 
and, in fact, present similar shape llBachallll979l . 

Despite the fact that the early descriptions of galaxy 
clustering properties were given in terms of the MF 
llGott & Turneiill977l) . most authors have focused on the 
shape of the mass function, which can be directly com- 
pared to the PS formalism. E yen when the observed quan- 
tity is the MF, some authors <jP>a,ctiall fc CenL ll9931 prefer 
to convert it into a mass function using a reliable M/L ra- 
tio. Nevertheless, one must consider all of the uncertainties 
introduced by the mass estimation, which are propagated 
to the mass function determination. These include errors 
in the internal velocity dispersion used for dynamical mass 
estimates, the large intrinsic scatter in the richness-mass 
relation, and errors in assuming dyn amical equilib r ium f or 
all clusters when using X-ray data ijGirardi et all fl99^l . 

The main problem which must be overcome when de- 
termining the MF is the production of a statistically signif- 
icant and unbiased catalogue of groups and clusters cov- 
ering a large enough area of the sky and encompassing 
cosmic structures spanning a wide range of richness, from 
very low multiplicity structures such as galaxy triplets, up 
to very rich clusters with several hundred members. 

In the past, catalogues of groups and cluster s have 
been derived f r om ei t her 3 - D dat a (cf. iMaia et all Q.989; 
iRamella et ail Il989l l200ll l2002f) or from projected (2- 
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iMaternd. Il978t Ide Filippis et all l2000|) . All these cata- 
logues are derived from different data sets and with dif- 
ferent algorithms and are therefore affected by different 
biases favouring the detection of structures in a given 
richness range; biases induced by the topology of the 
data, by the limited size of the survey, by a mbiguities 
in the selection criteria, etc. I Sheet marl l|l985l) pioneered 
the field of automated cluster finding in optical surveys 
using peak-finding methods, which has been refin e d and 



modif ied in many l ater projects Maddox et^lJiJl^^ 
Dalton et all Il992t ITumsden et all Il992t iNichol et al. 
2001aHGal et alll2002l) . Base d on a model dependent ap 



proach Tp«3tjrafflet*'ai (199G) developed the matched fil- 
ter techn ique, which has been widely used, with several 
variants feawasaki et all Il998t ISchuecker fc Bohringerl 
1998 I lLobo et all 120001). in cluding the adaptive matched 
filter IjKepner et all 11999^1 . In addition, the availabil- 
ity of multiband high accuracy CCD data, allowed the 
implementation of several cluster -finding methods based 
on the use of galaxy colours llCladders fc Yea . 12000: 
ICoto et all iH INichol et all 12001 ht lAudreorl M$ 
An independent approach relied o n the Voronoi tessel- 
lation technique as a peak finder ijRamella et all [2001: 
iKim et all UpOO) 

and a modified version, taking i nto ac- 
count colours was implemented bv lKim et alJ l|2002h . More 
recently, other, more advanc ed pattern recognitio n tools 
such as Bayesian clustering llMurtagh et all l2002h . max- 
imum likelihood ijCocco fc ScaramellaL Il999|) . and neural 
networks (Frattale Mascioli, Priv. Comm.) have been in- 
troduced. 

Much less work has been done to detect poorer 
structures such as loose groups; two principal methods 
(and their suc c essive elaborations) have been adopted. 
iTurner fc Gottl lll97fih presented the first tentative ob- 
jective identification of groups as enhancements above a 
reliable threshold in the projected gal axy distribution. 
The " Friends Of Friends" algorithm of iHuchra fc Gellerl 
(1982) generates a measure of correlation among galax- 
ies and their neighbours, basing on their separation in 
the full 3-D space. A noticeable exception to the lack of 
low-richness catalogs has bee n the detection of compact 
groups, where several teams ( de Carvalho fc Diorgovskii 



H995Hlovino et al] . ll999ll2002|) have proposed different ap- 
proaches to their detection. For the determination of the 
MF, it is important to note that its derivation from the 
above cited catalogues is hindered by the fact that all of 
the above algorithms are optimised for the detection of ei- 
ther groups or clusters, and no systematic work has been 
done in matching their outcomes in the transition region 
between structures of low and high richness. 

Here, we attempt the derivation of an accurate 
MF, starting from the galaxy catalogues extracted from 
DPOSS material. 

The paper is structured as follows. In Sect. |21 we 
shortly summarise the prop erties of the Digitized Palomar 
Sky Survey (DPO SS) data llDiorgovski et all ll998L fl999l 
Rei d et a,ll ll99l used to derive the multiplicity function 
described in Sect. El In Sect. OH we describe the algorithms 



used to detect groups (Sect. l3~T)l and clusters (Sect. l3~2*|) . 
while in Sect. 0] we discuss the simulations performed in 
order to evaluate the accuracy of the method, expressed 
in terms of completeness and fraction of spurious detec- 
tions, and to evaluate the possible existence of systematic 
errors in the ranges of overlapping richness for the group 
and cluster finding procedures. Finally, in Sect. [SI we draw 
our conclusions. Through this paper we assume Hq = 100 
km s~ 1 Mpc~ 1 . 



2. The data 

The data used in this paper were extracted from 
the DPOSS photo g raphi c plates l|Diorgovski et a 
Il999t iReid et all Il99lj) using the SKICAT package 
l|Weir et all |]J95a) which provides photometric, mor- 
phological and astrometric data for each detected ob- 
ject. SKICAT also provides a classification (Star/ Galaxy) 
based on a classification tree ijWeir et all Il995b[) . 
In DPOSS, the three photometric bands (J, F and 
N) are individually calibrated to the Gunn system 
llThuan fc GunnLll97fitlWade et allll979l) by means of ac- 
curate CCD photometry of objects of intermediate lumi- 
nosity, (to take into account the nonlinear response of the 
plates), with preferential targetting of galaxies. From the 
DPOSS data covering the selected regions, we extract, for 
each individual object: RA, Dec, total magnitude which 
best approximates the asymptotic magnitudes and the ob- 
ject classification. 

DPOSS individual plate catalogues must be cleaned 
of spurious objects and artifacts (such as multiple de- 
tections coming from extended patchy objects, halos of 
bright stars, satellite tracks, etc.). In order to do so, we 
mask plate regions occupied by bright, extended and sat- 
urated objects which locally make object detection ex- 
tremely unreliable. Subsequently, we matched catalogues 
obtained in each of the three photometric bands, by using 
the plate astrometric solution and by matching each ob- 
ject in one filter with the nearest objects i n the two other 
filter s (with a tolerance box of 7 arcsec, see lPaolillo et all 
2001). Due to the different S/N ratios in the three bands, 
many objects had discordant star/galaxy classifications in 
catalogues obtained in the different bands. The number 
of such objects obviously increases at faint magnitudes (it 
needs to be stressed, however, that this problem is greatly 
reduced whe n a new training set fo r the classification is 
adopted, see lOdewahn et alT l|2002') for details). In order 
to exclude from our final catalogues the smallest number 
of true galaxies, we discard only the objects classified as 
stars in all three filters. Final catalogues were thus ob- 
tained for 10 DPOSS plates (see Tab. covering a total 
area of ~ 300 sq. deg. spread at high galactic latitude 
(||6|| > 30 deg) (see Fig.^l, in order to reduce cosmic vari- 
ance. Details on the photome tric calibration of these par - 
ticular plates can be found in lPaolillo etHI ll200ll Eool . 
We note that these calibrations are no t the same as the 
general DPOSS calibrations described in lOal et alJ l|2f 
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Table 1. List of DPOSS plates, from whic h we extracted 
our ca talogues. Notes: (1) c alibration from JPaolillo et alJ 
l|200l|k (2) calibration from IPaolillo et all (|2003h . 
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Our catalogue of galaxies is limited in magnitude down to 
the Gunn r = 20.5ma<?. 

3. Detection of galaxy overdensities (groups and 
clusters) 

Making an arbitrary choice, we use the term "groups" to 
denote those galaxy aggregates which consist of less than 
20 objects, and "clusters" for all richer structures. This 
definition is comparable to that of Abell (1958), but in 
our case we set an implicit threshold on the magnitude 
difference between the brightest and the faintest objects 
in the same structure, given by the limiting magnitude. 

3.1. The procedure for groups 

In order to detect galaxy associations of low richness 
{N bj < 20), we have implemented a modified version 
of van Alba da's algorithm originally developed fo r binary 
systems (see lOosterlool Il989l ISoares et allTl995h . 

Taking into account only the position and the apparent 
magnitude for each galaxy in our catalog, we first search 
for the nearest neighbour in a given magnitude range, 
and then estimate the probability that the two objects 
are physically related. 

For the fore/background galaxies, the projected distri- 
bution is assumed to be Poissonian and the probability 
that the angular separation between a given galaxy and 
its nearest neighbour falls in the range 8 and 8 + d8, is: 



Pi(0)<20 = exp [-n8 2 p] 2n8pd8 



(1) 



where p is the surface density of background galaxies in the 
immediate neighbourhood. In order to combine the angu- 
lar separations of different pairs into a single distribution, 
the quantity x is defined as the ratio between the observed 
value of the distance {8\ ) to the nearest neighbour and the 
expected theoretical mean value < 8\ > given by Eq. ^ 



(2) 




Fig. 2. Upper panel: comparison of the Poissonian distri- 
bution (solid line) and the distribution of x's (histogram) 
in a simulation with galaxies randomly distributed in the 
sky. Lower panel, as upper panel, but for actual obser- 
vations: some of the nearest neighbours are physically 
linked (related) to the groups and produce an excess of 
ncibourghs at small x. 



The resulting frequency distribution of x: 

,•2 " 



Pi{x)dx = exp 



TTX 



-xdx 



(3) 



is then independent of the background density p. 



1yff> 



The shape of the observed distribution, pq{x), an d 
the Poisson distribution pi(x), for large x, are expected 
to be similar. If an excess is found in the observed 
distribution relative to the Poissonian expectation for 
small x (see Fig. lower panel), it is likely due to 
physical companions, which will tend to cluster at smaller 
distances than random projections. 

Normalising the observed distribution to the Poisson 
distribution, we can use the excess Pq(x) — Pi(x), observed 
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Fig. 1. Stereographic projection of a transequatorial sky region (ranging from S = —5 up to 20 degs and from R.A. = 
up to 80 degs) containing the ten selected DPOSS fields listed in Tab.^ In order to trace the global investigated area, 
all the detected clusters (without any cut in magnitude) are represented. 



at small x, to define the probability p that two galaxies, 
located at a certain x, are physically associated: 



p=l 



Pi 0*0 
Po(x) 



(4) 



In this formalism, all galaxy pairs having p higher 
than a given threshold value are considered to be physical 
companions. 

Iteration of the above procedure allows us to estimate 
the probability that other companions of higher order (up 
to N ~ 20) are physically related to the first object by 
comparing the observed distributions of higher order to 
the expected Poissonian distributions (normalised to the 
local density) for the second, third, etc. nearest neighbours 
(P2(x), P3(x), etc.). 

Groups are then identified by associating all galaxies 
having probability p higher than a given threshold value. 
Groups sharing one or more companions are finally merged 
into one single system. The total number of objects defines 
our richness for the groups. 

To compute the quantity x for every pair of galaxies, 
it is necessary to have an accurate estimate of the local 
galaxy density background p. To derive p, for each galaxy 
and within each magnitude interval, one first determines 
the distance to the i-th nearest neighbour 9i. The rela- 
tion between 9i and p is given by the probability that the 
distance to the i-th nearest neighbour lies between 9 and 
9 + d9: 



P t (9)d9 = exp [-tt9 2 P ] 



(i-1)! 



The mean expected value of 9i is: 

(e)= r(z + o.5) 

[(»-l)!v^rp)] 



(6) 



-2%6pdd 



(5) 



The higher the chosen value of i (i.e. for large distances 
to the i-th neighbours), the lower the probability of 
being affected by possible physical companions, which 
would lead to an overestimate of the local background. 
Furthermore the width of the distribution of the ratio 
between 9i and its mean value (< 9i >) decreases with 
increasing i. Thus if i is large enough, it is possible to 
obtain an accurate estimate of p from Eq. [S] by replacing 
the expected value < 9i > with the observed one 9i. On 
the other hand, i must not be too large, otherwise too 
much of the small scale clustering would be lost, and a 
large area of the plate will be affected by border effects 
(distant companions of galaxies located near the border 
of the plate will not follow a Poissonian statistics and will 
be preferentially located towards the center of the plate) . 
The choice of the value of i is therefore a compromise 
that has to be made by taking into account all ofthc 
above factors. 



3.2. The procedure of cluster identification 

Candidate clusters were ident if ied following a procedure 
similar to that of I Sheet man! <|1985l) . The catalogue of 
galaxies is binned into equal-area square bins in the sky, 
generating a two dimensional map (density map) of the 
number density of galaxies. The bin size (1.2') is chosen 
such that the mean number of galaxies per bin is ~ 1, in 
order not to degrade the spatial resolution. The resulting 
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Fig. 3. The smoothed two dimensional density map of the 
number density of galaxies for a field 5° x 5° centered at 
RA = lh and Dec = +15° (DPOSS plate n. 610). The 
smoothing has been performed by a Gaussian 2-D filter 
with width (~ 3') typical for a cluster core at redshift 
0.1 - 0.2. 



map (Fig. 0) exhibits irregular structures corresponding to 
the presence of overdensities emerging above the intrinsic 
fluctuations of the background distribution. 

The large fluctuations existing in the distribution of 
background galaxies are due to the non-uniform back- 
ground galaxy distribution. Once the density map has 
been created, the analysis of these maps poses similar 
problems t o those of classica l phot ometry, so we use S- 
Extractor ijBertin &: Arnoutel Il996|) for the detection of 
areas showing enhanced signal. S-Extractor is run on the 
density map searching for objects with a minimum de- 
tection area of 4 pixels above a global threshold of 0.4 
times the Poissonian background noise estimated from 
each plate using a background map. The evaluation of 
such background is a crucial step, strongly affecting the 
final richness estimate. The use of S-Extractor poses sev- 
eral problems (which cannot be trivially solved) since it 
is optimised to work on images with Gaussian statistics, 
while in density maps there are too few objects per bin, 
and they are distributed according to Poissonian statis- 
tics, thus making the background determination provided 
by S-Extractor unreliable. To circumvent this problem, 
we were forced to derive the background map in an alter- 
native way. We first divide the original density map into 
sub-images of ~ 1° x 1°, and then compute the Poissonian 
mean in each box, subsequently performing a fit with a 2- 
dimcnsional polynomial function of first order. We found 
a mean background density of 1640 per sq.deg. with a a of 
148 galaxies. In this way, we remove those spatial frequen- 
cies higher (the clusters) than the mesh scale length. At 



the estimated typical redshift in our sample (z = 0.1 — 0.2) 
this scale corresponds to a linear dimension of 9-15 Mpc. 
This map was then subtracted from the global frame be- 
fore running the detection procedure. 

The resulting density map was then smoothed in the 
detection step using S-Extractor with a Gaussian 2-D filter 
in order to match the cluster density profiles and, since we 
are searching structures with almost a Gaussian core, the 
filter width was chosen depending on the expected average 
apparent size for the cores (~ 250 kpc) of clusters in the 
redshift range [z = 0.1—0.2) probed by our data. We stress 
that the choice of the otimal parameters strongly depends 
on the characteristics of the specific data sets and needs 
to be tuned on the simulations reproducing the behaviour 
of true catalogues. 

The extracted parameters characterizing the detected 
overdensities are the density centroid in absolute equa- 
torial coordinates (J2000), the isophotal area above the 
threshold, the S/N ratio of detection, and the number of 
objects inside the isophotal area, which we use to derive 
(after the background correction) our richness parameter 
for the clusters. 



4. Outlines of the simulation 

In order to test the limits of our group and cluster 
detection procedures, we performed simulations over a 
region having the same area and the galaxy counts as 
one POSS-II plate. In this way we could estimate the 
shortcomings of our procedure, such as the percentage of 
spurious detections and the percentage of lost objects; 
at the same time this helped in the fine tuning of the 
parameters of the detecting algorithms. 



4.1. Simulation of the galaxy background 

First we simulated the galaxy background assuming a 
uniform galaxy distribution. The number of simulated 
background galaxies is the average number of galaxies 
present in the DPOSS plates (approx 50000 after exclud- 
ing all the galaxies fainter than the limiting magnitude). 
To each background galaxy, a sky position, randomly 
extracted within the plate limits, and an apparent 
magnitude, distributed according to the observed galaxy 
counts, were assigned. 



4.2. Simulation of galaxy groups 

The number of groups to b e simulated was extract ed from 
the multiplicity function of iTurner & Ootil l(l976h . 
We began by placing the principal galaxy of each group 
at random positions inside each field. Then, to each 
principal galaxy we assign an absolute magnitude and 
a redshift. Absolute magnitudes were extracted from 
a Schechter function with M* = —19.80 and a = 
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Fig. 4. MF of simulated (filled circles) and detected 
(empty squares) groups. On the horizontal axis there is 
the number of galaxies in each group, that is the richness. 



-1.25 jRamella et all [1999). while the redshifts were as- 
signed from the galaxy distr ibution observed in the Las 
Campanas redshift survey IjShectman et al 
each principal galaxy we then associate a number of sec- 
ondary galaxies matching the multiplicity function men- 
tioned above, each of these galaxies having the redshift of 
the corresponding principal galaxy. Taking into account 
the estimates provided in the literature, each simulated 
group was given a maximum standard dimension depend- 
ing on its richness: a maximum radius of 0.26 Mpc for 
groups with N tj < 18 members, while a maximum radius 
of 0.55Mpc is used for groups with N bj > 18 members. 
All the secondary galaxies belonging to a group were then 
distributed inside the group volume, and each assigned an 
absolute magnitude generated from the same Schechter 
function as the brightest galaxies in the group. Finally, ab- 
solute magnitudes were re-transformed to apparent mag- 
nitudes by taking into acco unt the cluster distanc e and the 
average fc-corrections from lFukugita et alJ 1 199.4 ). The de- 
tection algorithm was then applied to the simulated plates 
in order to fine tune the algorithm parameters (threshold 
value of the probability p and choice of the i — th near- 
est neighbour to compute the background galaxy density). 
The results of the simulations may be summarised as fol- 
lows: the group detection algorithm loses 28% of the sim- 
ulated groups and produces 43% spurious detections. 
Fig.0]shows that, in spite of the high contamination level, 
the MF shape is statistically preserved: the simulated MF 
(filled circles) and the detected MF (empty squares) differ 
on average by a vertical offset, which we take into account 
to correct the final group MF. 



In Fig. El (left panel), we show, as an example, the 
outcome of one typical simulation. The centers of the 
simulated (dark dots) and detected (empty triangles) 
groups are plotted; a circle is drawn when the two match. 



4.3. Simulation of clusters 

Cluster simulations were performed with same assump- 
tions used for the groups, with some crucial differences. 
The number of simulated clusters of a given richness (rang- 
ing from 2 to 200 galaxies) in an area of 37.59 squared de- 
grees (approximately the area of one DPOSS plate) was 
determined from a preliminary analysis performed on 10 
DPOSS plates. In a second step, a power law multiplicity 
function was used, with the slope taken from the prelimi- 
nary multiplicity function. In this way we tried to take into 
account the total number of low richness objects, which 
could not be measured from our preliminary analysis. 

The absolute magnitudes of the principal galaxies 
were extracted fro m a Gaussian distribu tion centered on 
-22.99 ± 1.0 mag Schneider et allll983ft . while those for 
the secondary galaxies were extrac ted from the luminos- 
ity function of IPaolillo et al.l l|200l|) . To take into account 
the richness dependence of the cluster dimensions, we ar- 
bitrarily adopted a core radius (<r of the Gaussian profile) 
of 0.5 Mpc for clusters with < 30 members, while a core 
radius of 1.0 Mpc was used for clusters with > 30 mem- 
bers. Although these values may appear somewhat high, 
the adoption of smaller values for the core radius would 
only make the detection easier and therefore the whole 
procedure more reliable. As with the groups, the detec- 
tion algorithm was applied to a large number of simulated 
plates to test the algorithm performance as a function of 
the properties of the objects to be detected. 

In Fig. El we plot, for a typical simulated plate, the 
assigned richness vs. the assigned core radius of the sim- 
ulated clusters (open circles) and mark with a cross the 
clusters retrieved by the algorithm. Clusters with a very 
shallow profile or which are poor are preferentially lost. 

The dependence of the algorithm efficiency on the rich- 
ness is shown in Fig.0 where we plot the number of sim- 
ulated (continuous line) and retrieved (dash shaded area) 
clusters in the typical plate area. All but two of the clus- 
ters having N j,j > 35 are retrieved. In the range of rich- 
ness 25 < N bj < 35, 80% of the clusters are retrieved. 
Considering that a cluster belonging to the Abell richness 
class (30 — 49 members in a range of two magnitudes) has 
N bj > 30 (N bj includes the cluster galaxies in a range of 
at least four magnitudes), we are complete up to z = 0.2 
at least for all the Abell richness classes. Fig.|J|also shows 
that spurious detections (dot shaded area) are absent in 
the richness range where the algorithm works with the 
highest efficiency, and occur only in the range where the 
group finder is to be used. 

As already mentioned, the estimate of clusters richness 
is given by the number of objects within the detection 
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Fig. 5. Simulated (dots) and detected (triangles) structures. Left: groups. Right: clusters. Circles highlight simulated 
groups/clusters which have been detected. 
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Fig. 6. The richness vs. the core radius of the simulated 
clusters (open circles). The crosses mark each retrieved 
cluster. 
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Fig. 7. Cluster number in the one plate area is plotted 
as a function of estimated richness. The continuous his- 
togram represents the number of clusters given in input 
to the simulation; the dash shaded histogram represents 
the retrieved clusters and the dot shaded histogram the 
spurious detections. The richness bin grows exponentially 
as 2™/ 2 (see Sect. O . 



isophote (isodensity counts). We wish to stress that 
this definition of richness depends on the redshift of the 
detected structure. The quality of the richness estimate 
has been tested using our simulations. In Fig. |H1 the 
points follows bisector of the diagram (a bit shifted 
towards the upper half part of the plot), with a scatter 
in richness of ~ 10 galaxies (which is consistent with 
the background fluctuations). The small shift indicates 
an understimation of the retrieved richnesses. We are 



comparing the number of the galaxies put in a synthetic 
circular aperture (the simulated) with the richness in the 
isodensity irregular countours, as it is measured in the 
real case: in this way some galaxies are missed. If we use 
circular apertures of the cluster size (which are known in 
the simulations but not in the actual observations), the 
shift disappears. Points in the lower right part of the plot 
are due to overlapping clusters, for which (in the absence 
of a deblending procedure) the richness will obviously be 




Fig. 8. Richness of the simulated vs. the detected clusters. 
The errors are inversely proportional to the signal to noise 
ratio for the detection. 

overestimated. 



Fig. 9. The multiplicity function for galaxy structures 
ranging from small groups (filled triangles) to rich clusters 
(filled circles). We remove clusters in the richness range 
where detection efficiency is low. 



5. The conjoined groups/clusters multiplicity 
function 

Fig. EH summarizes our main results. We plot the MF, de- 
fined as the number of groups or clusters per unit area and 
per unit of estimated richness (the groups/clusters rich- 
ness is defined respectively in Sect. IH. II and Sect. l8.2|l . For 
the clusters, the bin grows exponentially as 2"/ 2 , in order 
to keep the S/N ratio almost constant along the richness 
axis. For the groups, the bin was instead set equal to 1. 
In order to exclude the structures detected in the redshift 
range where our magnitude limited catalogue is incom- 
plete, only clusters and groups where the brightest galaxy 
has m < 16.5 (in Gunn r) were selected. Assuming that 
brightest galaxies may be used as standard candles, our 
selection in magnitude implies that z < 0.2. 

The procedures described above were applied sepa- 
rately for groups and clusters, obtaining two different mul- 
tiplicity functions (marked with different symbols in Fig. 
|5J). These MFs appear to define a common relation, with- 
out the need for any offsets or normalisations. We empha- 
size that a minor correction for completeness was applied 
only to the last point of the clusters MF. To correct the 
group's MF for contamination by spurious detections (see 
Sect. IHHl , a global shift derived from the simulations was 
also applied. Only the Poissonian statistical fluctuations 
have been taken into account in the error estimate. For 
the high richness clusters, the error on the richness esti- 
mate is negligible with respect to the bin width. The error 



becomes relevant only in the same richness range where 
incompleteness is also significant. 

In Fig.^]]we compare our results with a MF extracted 
by us from the USGC catalog of groups I Ramella et all 
2002). We adopt the same representation scheme for the 
two data set. Normalisation to the same volume was ap- 
plied to the USGC groups, assuming a uniform distribu- 
tion of objects in redshift both for our sample and the 
USGC sample. It is important to note that the two cata- 
logs were derived in totally different ways. The USGC is 
generated using spectros copic redshifts by a p ercolation 
method implemented bv iRamella et alJ Jl997l) for group 
detection, which is designed to reduce the risk of false de- 
tections introduced by chance projections. 

The agreement between these two MFs (see Fig. HOf) . 
derived under totally different assumptions and using in- 
dependent data sets, is due to similar biases affecting the 
estimated richnesses for both samples. For low N bj struc- 
tures (groups) the similarity is apparent; in both cases, the 
methods count individual objects fulfilling the respective 
membership criteria but with secondary members having 
magnitudes falling within similar (i.e. ~ four magnitudes) 
ranges with respect to the primary galaxy. For the clus- 
ters, instead, the different depths sampled by the two data 
sets, when compared to the different limiting magnitudes 
of the samples themselves, indicates that both methods 
sample very similar intervals of the cluster's luminosity 
function. 
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Fig. 10. Overplot of the MF of USGC2 groups (empty 
circles) on the multiplicity functions obtained from the 
DPOSS data. 



6. Summary and discussion 

We have implemented two algorithms for the detection of 
galaxy associations, one for groups and one for clusters. 
The former is a modified version of van Alba da's proce- 
dure t o detect galaxy pairs, while the latter the lShectmanl 
l)l985|) approach, which uses a peak-finding procedure on 
a density map obtained from the galaxy catalogue. 

We evaluated the performance of these methods 
via extensive simulations, which show that the group 
algorithm is reliable up to richness 20, and the cluster 
algorithm is reliable at richnesses above 20 galaxies. 
The two algorithms were then applied to a ~ 300 square 
degree field extracted from DPOSS data (see Sect. O ■ The 
resulting MFs show a remarkable internal consistency 
from the two procedures which produce independent 
MFs for groups and clusters, matching with no need for 
normalisation. Additionally, the MF deri ved using our 
techni que on the 3-D based catalogues of iRamella et al.l 
l|2002h agrees with the MF derived from the projected 
DPOSS data. The final combined MF is well fit by a 
power-law of slope a — —2.08 ± 0.07. The correlation 
coefficient on the log-log scale is —0.98. 

The data set we used to determine the MF samples 
a volume [~ 300deg 2 , z < 0.2] whic h is sli ghtly smaller 
than that explored bv IBachall et all l|2002l) [~ AOOdeg 2 , 
z = 0.1 — 0.2]. The total number of d etected structures for 
N > 10 in the IBachall et alJ ((2002) and in our sample is 
respectively ~ 300 and ~ 370. In a forthcoming paper we 
shall analyze the cosmological implications of the derived 
MF. 
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